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Abstract 

We have calculated the asymptotics of Green's function of the differential equa- 
tion of nonlinear diffusion in the microscopic range with strong porosity fluctuations 
in the problem of flow through irregular packed beds for the arbitrary dimension of 
space and arbitrary porosity fluctuations covariance. 

PACS numbers:47.55 M, 11.10.G , 47.10. +g, 05.40. +j, 42.27. Gs 

1 Introduction 

Properties of flows in porous media are of great interest of modern industry and technol- 
ogy. In spite of practical importance and long-standing attention paid to this problem 
there is still no solidarity on model of the flow through the packed beds. 

One can point out an approach of a single particle diffusing with random dynamics in 
a Poissel velocity field, which goes back to the classical work of Taylor . Nevertheless, 
the experimental measurements are still described by the usual diffusion equation where 
a modified dispersion coefficient is used; see and references therein. The problem of 
these considerations is of the essentially gaussian solutions of the diffusion equation. The 
numerous experimental observations, || demonstrates clearly that the previous approach 
needs to be supplemented. 

One of the first attempts to overcome this limitation is proposed in || , where a Markov 
random process is used to model the diffusion of impurities in a turbulent stream. In [[| the 
stochastic field was introduced to model the porosity instead of the diffusion process itself. 
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Their result for the one- dimensional Fick law of diffusion for some marker transported by 
a carrier fluid through an irregularly but statistically homogeneous packed column leads 
to the appearance in the equation of the mean concentration of a time diffusion term. 
Together with a source term they are to be responsible for the departure of the solutions 
from gaussian. 

In this paper we keep on study the diffusion equation of the Fick-type in which the 
porosity e fluctuates strongly along the column axis. There are many methods of defining 
porosity; see || and references therein. Each of them replaces the complex network of 
voids with a single number that represents an average property. In the present work we 
use porosity as a point quantity in a porous medium and define an elementary volume to 
characterize a medium macroscopically following || and J7|. The size of the elementary 
volume around a point P we suppose to be smaller than the total medium so that it can 
represent a fluid flow at P. We suppose also that there is enough pores to allow statistical 
averaging. Since porosity varies, the maximum length is the characteristic length £ that 
indicates the rate of change of porosity, and the minimum length w is the pore size. 

The porosity is directly related to the size of the pores relative to the matrix. Consid- 
ering the problem of packed column, we denote by x the axial coordinate of the column, 
< x < L, then the fraction of the column that is voids defines the averaging value of 
£, < e < 1. If P be a point inside a porous medium on the column axis, see Fig.l., 
surrounded by a layer of volume V(/j), one can define a ratio 

_ Vy{k) 

£l V(k) 

in which V v (li) is the volume of void space within volume V(k). 

Consider a sequence of values of k such that li > 1% > l 3 > For large values of 

V(li) the ratio £j may change gradually as V(lj) gets smaller. As I falls below a certain 
value £ there will be large fluctuations in porosity when V(l) is approaching the dimension 
of the pores (this V(^) may be taken as the elementary volume). Below this value the 
average value of porosity has no meaning, then it is natural to consider the porosity as a 
function of x, where i be a distance apart from the point P. 

Flow phenomena for x > £ is well investigated in the porous medium theory; see for 
example for a review. Below £, for w < x < £ flow properties are subject to strong 
microscopic effects and usually left beyond the investigations due to an extremely large 
fluctuations in porosity when x is approaching the dimension of the pores. However, 
namely this interval of scales is of great interest for the variety of problems in chemical 
engineering and environmental studies. 

We continue to study the microscopic properties of flows through a porous media 
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with strongly fluctuative porosity using the quantum field theory renormalization group 
method (RG) and 4 — e expansion ||. The technique which is developed in what following 
is similar to one of |§ devoted to waves propagation in a randomly inhomogeneous medium 
with strongly developed fluctuations. 

In Sec. 2, the problem is stated and formulated in the quantum-field formalism. In 
Sec. 3, the Green's function for the equation of nonlinear diffusion is constructed for the 
case of gaussian distributed fluctuative porosity. We have introduced the effective action 
functional which is equivalent to the problem discussed under the certain conditions in 
Subsec. 3.1. Investigating the properties of the theory, we have justified the result of 
||] on the inclusion of a "time-diffusive" term into the equation of nonlinear diffusion in 
case of fluctuating porosity for arbitrary dimension of space and for arbitrary stochastic 
process originating the porosity fluctuations. 

In Subsec. 3.3 and Subsec. 3.4, the procedure of ultraviolet renormalization of the 
theory with the "time-diffusive" term included is considered and critical dimensions in- 
cluding anomalies are calculated. The results obtained for static correlation functions 
meet the well-known empirical laws of Richardson and Kolmogorov. It is shown that for 
the parameters taking values within the physical domain (positive viscosity, etc.), the 
renormalization-group equations have a manifold of infrared (IR) attractive fixed points. 
This implies the presence of a universal scaling regime in the IR domain of the micro- 
scopical region of fluctuative porosity for some range of values of physical parameters. 

The main result of the present paper is the detailed description of the effect of strongly 
fluctuating porosity on the form of static ( i.e., of equal time) as well as of dynamic Green's 
functions of the differential equation of nonlinear diffusion. 

Similar to the problems of nonlinear diffusion and of chemically active scalar admix- 
ture, the problem of the flow through the irregular packed beds with strongly fluctuative 
porosity, in the general case, is of an infinite number of coupling constants, i.e., no 
one of the statistical momenta of higher order can be omitted from the consideration. 
However, in contrast with the problems mentioned any of these terms does not contain 
UV- divergences, so that the UV renormalization procedure as well as the critical behavior 
of the quantities are unsensitive to certain hypotheses on the porosity function as well as 
the certain statistical distribution for porosity in the microscopical range. 

It turns out to be that fluctuative porosity does not affect on the universal properties 
of the static Green's functions, i.e., it does not change the indices of empirical laws 
of Richardson and Kolmogorov. But porosity changes the amplitude factors of static 
correlation, i.e., the Kolmogorov constant and the effective diffusivity coefficient. 

Fluctuative porosity does affect the dynamic Green's functions providing the damping 
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retarded time-spectrum. 

Qualitatively, for the concentration one has an outgoing damping traveling wave type 
solution which decay rate as well as an amplitude depend on porosity. 

2 Formulation of the Problem 

We shall consider the diffusion process of mixture in a randomly inhomogeneous packed 
beds with strongly fluctuative porosity in the microscopic range w < |x| <C £ for the 
arbitrary value of space dimension d. 

The problem is posed as follows: To find the retarded Green's function (propagator), 
averaged over the imposed statistics of fluctuative porosity e(x) and flow velocity v(x, t), 
of the inhomogeneous diffusion equation for the concentration «(x, t) of some marker in 
the available fraction of the volume: 

£(x)<9 t u(x, t) + (vd)u(x, t) + z/Au(x, t) = (1) 

in which v is the Fick's diffusion coefficient, A is the Laplace operator, v(x, t) is the 
velocity of fluid flow. The porosity e(x) = £ + £( x ) is the sum of constant and random 
components. From the physical point of view the value of porosity belongs to the interval 
s E [0, 1] , so that, generally speaking, e cannot be taken as a gaussian distributed function. 
To model porosity one can use any function which values lie within the interval [0, 1] ; for 
a qualitative determinacy in what following we model it by the function 

e(x) = (l + sinc-V(x))/2, (2) 

where c = is a parameter with the dimensionality of velocity which characterizes the 
velocity of the diffusion process in the scale £, and <£>(x) is a gaussian distributed phase 
function (of the same dimensionality) with zero mean value and static correlation function 
of the Ornstein-Zernike type, |l0j 

D v = MxMy)> = Vo^jAa / dk(k 2 + f exp ik(x - y), (3) 

where \l ~ l/£ is a mass parameter. We consider all parameters in (|3|) as bare (non- 
renormalized ) parameters and supply them with the index "0". For future convenience 
in (^) we have introduced the amplitude factor by such a way that it meets the proper 
dimensionality of the correlation. The coupling constant t]q plays the role of an expan- 
sion parameter in the perturbation theory such that the relation t/qUq = A 2S defines the 
characteristic maximum momentum scale in the theory, h.~l/w. 
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The parameter of the regular expansion of correlation functions 25 is the deviation of 
the space dimensionality from 2. We demonstrate below that under the certain physical 
assumptions the real value of parameter 5 is S r = 1. 

The dragging velocity v along with the fluctuating component of porosity e prescribe 
a certain dimensionless parameter in the problem which is analogous to the Reynolds 
number, Re £ = evw/u, where e is characteristic porosity fluctuations in the scale of w, 
and v is a characteristic fluctuative flow velocity in a pore. A statistical steady state is 
expected when Re e ^> 1. In terms of the energy dissipation of fluid flow the parameter Re e 
defines a dissipation length Id by the relation Re^J 3 ~ IdW^ 1 . Starting from some values 
Re £ , the dissipation length falls into the range of microscopic effects, w <C r <C £. It leads 
to the formation of inertial range of scales, Id <C £, which is substantially analogous to 
the Kolmogorov's dissipation free range |16| . 



The basic results of the phenomenological Kolmogorov-Obukhov theory JT^| can be 
generalized directly to the problem in question: Assume that in the inertial range of 
wave numbers (l/£ <C k <C 1/Z^) correlation functions of (p depend not only on the mean 
energy dissipation rate W of fluid flow but also on an additional argument, the velocity 
of diffusion process, c. Then, assuming that the correlation functions do not depend on 
viscosity, for an equal-time pair correlation function of ip, one has 

D{k) ~ W^k-^TiW^k-^c- 2 ), (4) 

where T is a scaling function of dimensionless arguments. By the way, ([|) predicts the 
Kolmogorov's scaling (A v = —1/3) for the stochastic phase field (p. This phenomeno- 
logical result is taken into account in the model correlation (|3]) by choosing 5 r = 1 as 



a real value for 8. In accordance with m) and |]I8 |, one obtains the scaling dimensions 
for quantities in the amplitude factor of © A c = —1/3 and A u = —2/3; for 5 = 1 (||) 
provides exactly the Kolmogorovian scaling dimension for the correlation (<p(p). 

Similarly, to describe the effect of eddy diffusion which obviously plays an essential 
role in the microscopic diffusion process we consider the dragging velocity to be a d- 
dimensional gaussian vector field v(x, t) of zero mean value and of covariance: 

D v = ( Vi (x, t) Vj (y, t')) = 9o^ t~J J dkP i:j (k)(k 2 + m 2 )-$-iexpik(x-y), (5) 

where Pij(k) = dij — kikj/k 2 is the transverse projector (the fluid flow is considered as 
incompressible), go is a relevant coupling constant ( another expansion parameter in the 
perturbation theory), 1/m is a turbulent integral scale A' which can be chosen as A' = A; 
< e < 2 is another parameter of regular expansion of Green's functions which real 
("Kolmogorov's") value is e r = 2/3. Note, that the velocity correlation (0) as specified 
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above is inherent for numerous simplified models of turbulence introduced by Obukhov, 
II | and Kraichnan, [12], see jnj] and references therein. 



The object of our interest is G = ({G(<p, v))), where G((p, v) is the Green's function of 
the differential equation (p]) for p and v fixed and the doubled angular brackets denote 
the operations of functional averaging over the known statistics of p and v. 

The stochastic phase factor sinc~V(x) can be expanded in powers of <p, 



sine V( x ) = X! 



(c-V) 2 "- 1 ! 



n =i 2n - 1! 

as well as the averaged product of two phase factors is expressed by a series of various 
statistical momenta, 

» ^ ((c-V) an - 1 (x)(c- 1 y) am - 1 ( y )) 
(sine ^(x)smc <P&)) = ^ 2 n-l!2m-l! " 

In order to sidestep consideration of the onset of the problem we compute first the Green's 
function of the simplified equation with no terms proportional to (p (x), k > 1. Contribu- 
tions of the higher statistical momenta then are to be determined by various composite 
operators Fk = (c~ k ip k (x)) which we study in the forthcoming sections. 

We conclude this section by an explanation of the physical meaning of propagator G. 
The Fourier transformed propagator G(k), k = |k| is to be computed in the framework of 
diagrammatic technique of Feynman. Assuming that initially the marker was concentrated 
in the source of frame, J(x, t) = 5(x)/(i), one obtains the relevant concentration profile 
of the traveling wave type, u(r, t) = u(r)f(t — re/v) by the expression 



u[r) 



(2tt) - d J dxG(k) exp zkx. (6) 



Another practically important problem is to define a concentration distribution profile in 
a porous half-space. The flat boundary of the medium plays the role of a source term in 
the right hand side of the equation (|I|). In case of the porous half-space x\ = r > 0, the 
concentration profile (u(r, t)) inside the porous medium is determined by the equation ([I]) 
with a source term of the form J(x, t) = 5(x\)f(t) and yielded by the same formula (|6]) 
with d — 1. 

Notice that (|6]) allows one to bring about the angular integration for arbitrary value 
of space dimension d, 

POO 

u (r) = (27i)- d/2 r 1 - d/2 / dk G{k) k d/2 J d/2 ^{kr). (7) 



where J u is the Bessel function. Qualitatively, the solution of (0) is always an outgoing 
damping traveling wave of concentration which decay rate is of essentially interest for 
practical applications. We shall discuss the calculation of this rate below in Section 4. 
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3 Solutions for the Case of Gaussian Distributed e 

To illustrate our approach we start our consideration with unphysical but the simplest case 
of gaussian distributed fluctuative component of porosity e. In the forthcoming sections 
we recover the results for the general case. 

3.1 Diagram Technique and Effective Action Functional 

In the present section we develop the diagram technique relevant to the problem discussed. 
The bare propagator Lq 1 (the Green's function of the differential equation ([!]) with no 
nonlinear terms proportional to <p and v) is a retarded function, 

(^(k, t)^(-k, 0)o = L \k, t-t') = 6(t - Oc"^^, (8) 

and can be naturally pictured as an oriented line in diagrams. Diagrammatic expression 
for the propagator G(ip, v) is 



G{<p,v) = 




in which the dashed tail corresponds to the field v and the saw-type tail is associated 
with the field (p. In what following functional averaging over (p and v the relevant tails are 
connecting to each other, producing the dashed lines for the velocity correlation function 
(|) and the saw-type lines for the stochastic phase correlation function @. Up to one- 
loop order the diagram series for the averaged propagator G is forthcoming. The diagrams 
are similar to the diagrams of ordinary quantum field theory in which some "physical" 
complex valued field ^(x, t), ^(x, t) supplied by the bare propagator Lq 1 interacts with 
a couple of stochastic fields y(x) and v(x, t). 
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G(v>,v) = 




The relevant action functional has the form 

1 „ n 1 



S(<p,v,il>) = |dxdt --ipD~\ - 2 vD v lv + ^ ( L oV> + ^ + (v0)VO 



(9) 



C.C. (10) 



where C.C. denotes a complex conjugate part, -D^, and _D„ are the correlation functions (|3]) 
and (|) consequently. Averaging the product ^(x, t)if)'(y, t) over the complete set of fields 
{ip, v, for the statistical weight exp 5'(<^, v, one computes the dressed propagator 
(the correlation function G) of the "physical" fields, 

G(x, y) = «V(x, t)^(y, t))) = C J DipDvT)ip ^(x, t)^(y, t) exp S(y>, v, # (11) 

in which the factor C is defined by the normalization condition ((1)) = 1. 

The correspondence between the original problem ([]]) and the field theory ( |IUD is not 
exact, since in the field propagator (]TTf) there are redundant diagrams containing closed 
loops of the ipip'-lxnes. These diagrams were contributed to the correlation functions 
((p(p) and (viVj) which we think of to be exact as defined in (|3|) and (0). In principle, 
these diagram are to be eliminated; one can simply say that the original problem is 
equivalent to the field theory ([[(]) except the diagrams containing the closed loops of lines 
of "physical" fields ifjifj^. This elimination does not affect the theory, since, first, we are 
not going to perform any functional variables transformations, which could mix up the 
redundant diagrams with those we need, and second, the renormalization group perfectly 
well permits such an elimination of a certain class of diagrams. 

We also note that all closed loops of ^^-lines with no time derivatives in a cycle are 
automatically equal to zero in the theory (|T0|), since it contains retarded and advanced 
functions. In particular, from the fact of elimination of all diagrams containing closed 
loops of oriented lines it follows that 

(i) correlation functions of any number n > 2 of the stochastic field (p are trivial, 

(p(xi) . . .y(x„)) = 0; 

(ii) correlation functions of any number n > 2 of the stochastic field v(x, t) are trivial, 



(v(x 1 ,t)...v(x n ,t)) = 0; 
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(iii) any mixed correlation functions of stochastic fields y(x) and v(x, t) are trivial, 

( ¥ >(xi)...v(x n ,t)> = 0. 

Furthermore, in the model (0), the odd multipoint correlation functions of scalar field 
vanish, while the even single-time functions satisfy linear partial differential equations. 
The solutions for the pair correlation function can be obtained explicitly in analogy with 
that of in the Kraichnan model of turbulence [ p^j . 

Finally, it is essential that the action (PH|) is invariant under the following field trans- 
formations: 

V>(x, t) -> ip(x + s, t + r), cp b (x, t) -> ip(x + s,t + r) - b(x), 
v a (x, t) -> v(x + s, t + r) - a(t), 

with two parameters: an arbitrary vector function of time a(£) decreasing at t — > — oo, 
s(t) = ftoo a(t')dt' and an arbitrary scalar function 6(x) fading out at |x| — > oo, r(x) = 
/^>|x-x'[>w b(x )dx'. 

The integral turbulent scale as well as the dissipative length are taken into account 
in the model ( |T0D by the parameters of infrared cut off m and ultraviolet cut off A. The 
Ward identities which express the invariance of the model ( |i0|) under the transformations 
(|l2l) guarantee that all UV singularities of diagrams are subtracted out in each order of 
perturbation theory, so that all correlation functions have finite limits in — > and do 
not depend on in the inertial range. In the region m — > there are strong infrared 
singularities (powers of m) in the perturbation theory. However, these singularities are 
not related to the dynamic interactions of eddies which form the spectra of the propagator 



G but to the kinematic effect of dragging of small eddies by the large ones [ 13[] . This effect 
can be eliminated completely by taking of the frame of reference which moves with an 
arbitrary speed of the large-scale eddies. 

Since these singularities do not contribute to the spectrum for G, we will not take 
them into consideration, assuming that the frame of reference moves with some arbitrary 
speed v along the direction of flow. By the way, the solutions for concentration is to be 
a traveling wave with the time argument t — re/v, as it was proven in H. 

3.2 Dimensional Counting. Action Functional with a "Time- 
diffusive" Term Included 



We start now with the dimensional analysis of model flTD|). Since the model fliTf ) is of two 
scales, one can introduce two independent canonical dimensions to each quantity F in the 
theory (the momentum dimension, dp, and the frequency dimension, d F ). Then, based 
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on d k F and dp, on can introduce a total canonical dimension dp = d k F + 2rf£. (in the free 
theory, ~ A). Assuming d\ = d^ = 1, one finds out the dimensions of quantities in 
(|l0l) (see Tab. 1). Nonrenormalized parameters in the Tab. 1 are supplied by the index 
"0". 

Superficial ultraviolet (UV) divergences, whose removal requires counterterms, can 
present only in those Green's functions T, which canonical dimension 

5 v = d+2-d 4> N 4> (13) 

is a nonnegative number, ||15|| . Here, is the number of corresponding fields entering 
the function T and <p = {if, v, tp}. 

In the model ( JlOj ) the derivative d at the vertex ipiyd)^ can be moved onto the field 
ip by virtue of the vector field v is transversal. Similarly, the time derivative dt at the 
vertex c~ l ip(pdtip^ can be moved onto the field ip by virtue of the field ip does not depend 
on time. It decreases effectively the real index of divergence: 5p = 5r — N^. 

From the dimensions in the Tab. 1, taking into account the auxiliary considerations 
(i)-(iii) of the previous Section as well as those of discussed above, one can find out that 
for any d superficial divergences can only exist in the 1-irreducible function (ipip ') for 
which 5r = 2, 8^ = 0. The corresponding counterterms must contain either two symbols 
d and is therefore reduced to i/jAi/;* or two time derivatives along with the squared c" 1 , 
vc~ 2 ipd 2 ip^ . Inclusion of these counterterms into the action functional (|1(]) cannot be 
reproduced by a simple multiplicative renormalization because of the latter counterterm 
is not formally present in the action functional fllPl) . 

To ensure the multiplicative renormalizability, we introduce an additional term z/ Cq 2 x 
ipd 2 ip^ into (fl0|) supposing cq to be a bare parameter; it leads to a renormalized action 
functional of the form: 

S R ((p, v, V) = / dxdt [-i^D- V - |vD v x v + ^ (d t ^j + Zxz/A^ + Z 2 v C - 2 dfy+ 

+ c- x <pdtil) + {vd)^)\ + C.C. 

(14) 

The inclusion of the new term proportional to c -2 corresponds to adding a term of the 
type uc~ 2 d^u('x,t) into (p. Notice that the inclusion of the z/c -2 <9 2 u(x, t)-term (a "time- 
diffusive" term) into the diffusion equation in the presence of axial depending porosity for 
the first time was proven in || for the case of particular stochastic processes originating 
the large scale porosity fluctuations. Our approach allows to justify this result of 0] 
and to generalize it for any stochastic law for a porous structure. The coefficient vc~ 2 is 
analogous to the time parameter r which was introduced in M to characterize time of the 
diffusion process measured in the time scale of the carrier fluid. 
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Taking into account the new term in the diffusion equation ([I]), one obtains instead of 
in the limit (k£) 2 > 1 

(V>(k, t)^(-k, f))o ~ c e- cklt '- tl /2uk. (15) 



The damping of correlations in (15) as c — > (£ — > oo) is an essential feature of the 
discussed approximation. 

However, we may investigate the initial theory with no such an addition, if we solve 
the RG equations obtained in extended model with the initial condition Cq 2 = 0, i.e., 
if we suppose that either the diffusion process goes up extremely fast or the correlation 
length £ is very short due to high dense packed bed. In both cases the domain of discussed 
microscopic phenomena happens to be small enough and the fluid flow behavior in the 
packed bed column may be well investigated within the framework of usual porous medium 
theory dwelling on the constant averaged porosity. 

In presence of the time-diffusive term the solutions (0) exhibit a decay in time which 
rate is determined by the singularity of propagator G, k 2 = % 2 , where \ 2 = v^ 2 /c 2 in 
Fourier representation is analogous to the squared relative refractive index in optics. In 
a simple pole approximation for G~ x = k 2 — \ 2 the Bessel function in ([?]) can be replaced 
by its asymptotics at large kr that yields the large distance asymptotics for u(r): 

u ( r ) = _L(_^_)(«*-D/2. 

2x lixir 

In this case Imx defines the extinction coefficient Im% ~ l/lo, where Iq is the damping 
length. 



3.3 Renormalization-Group Equations, Critical Scaling, and 
Empirical Laws of Kolmogorov and Richardson 

The UV-divergences (in our case poles in e and S in diagrams) of the extended model 
considered (with a "time-diffusive" term included) are removed by the multiplicative 
renormalization procedure. It amounts to the following: the initial action 5 , ($) is referred 
to as nonrenormalized, its parameters and coupling constants are referenced to as bare; 
these are considered as some functions (remaining to be determined) of new renormalized 
parameters and coupling constants. 

The renormalized action functional ( |i4|) is a function of renormalized coupling con- 
stants and parameters: 

c = cZ cl g = gM e Z g , 
7] = riM 2& Z^ v = vZ v , 
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where all renormalization constants Z a are the functions of two independent quantities 
Zi and Z 2 : 

Zy = Z 1 , Zg = Z 1 1 , 

7 _ v-5 7 2 7 _ 7 l/2 7 -l/2 K X ' ) 

~ 1 2' — ^1 ^2 ) 

which can be calculated within the framework of diagram technique. We choose the 
simplest form of subtraction scheme where the divergences are presented as the bare 
poles in regularization parameters e and 5 (so called minimal subtraction scheme (MS)); 
M is the renormalization mass parameter, g, rj, and v are renormalized analogues of the 
bare parameters go, r] , and u . Z a = Z a (g,r],c,e,S,d) are the renormalization constants. 
The relations (16) and (fl7|) result from the absence of renormalization for the nonlocal 



contributions of D v and D v in the action functional (|14|), so that g u = gM e u and 
t/qUq/cq 2 = M 2& r\v z I c~ 2 . No renormalization of the fields ip and v and "masses" are 
required, i.e. Z v = Z^ = 1 and M = M, m = m, Z m = Zm = 1- 

The only field tp requires renormalization, ip = i^rZ^, Z v = Z\^ 2 Z 2 1//2 , so that the 
renormalized Green's function W R meet the relation 

W R {g, V , c, u, M)Z»* = W{g , Vo , c , u Q ). (18) 

W R are UV-finite function ( they are finite in the limits e — > 0, 5 —>■ 0) for fixed parameters 
a. 

The RG equations are written for the functions W R which differ from the initial W 
only by normalization and, thus, can be used equally validly for critical scaling analysis. 
To derive these equation one can note that the requirement of eliminating singularities 
does not determine the functions eo = eo(e, e, 5), e = {g, rj, c, u}, uniquely because of the 
value of M is not fixed by any physical condition. Variation of M for fixed values of bare 
parameters eo leads to variations of e and Z e , Z v . Following the standard notation, we 
denote by Dm the differential operator M&m for fixed eo- Applying it on both sides of 
( fL8l) leads to the basic RG equation 

[D M + P 9 d g + (3 c d c + fadr, - lv D v + n vlv \ W R = 0, (19) 

where we have used D x = xd x for any parameters of the renormalization theory; for any 

7i = D M In Z t , (3 a = D M a, a = {g, r), c}. (20) 
These identities determine the (3— functions of the theory considered, 

Pg = -9 (e - 7i) , 

P v = -r ] (26-5 ll + 2 l2 ), (21) 
Pc = -f (7i - 72) • 
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and the anomalous dimensionalities 7$. One calculates the renormalization constants Z\ 
and Z2 from the diagrams of perturbation theory (for example, up to one-loop order from 
the diagrams depicted in (§)) and then 7, , j3 a functions. All 7- and (3— functions are 
constructed as series in g, r], c, and the functions 7^ depend on neither e nor S. 

Actually, the only reason for calculation of renormalization constants in the theory 
considered is to prove the existence of infrared (IR)-attractive fixed points {g*,?7*,c*} of 
the RG equations ( ]T5| ) in the space of renormalized charges such that /3 a (g*, V*> c *) = 
and u a Q = <9 a /^(g*, r/*, C*) > 0. The values of 7^(7*, 77*, C*) then can be found from (|2l|) 
exactly apart from diagram calculations, which are, in fact, pretty standard. 

There is a line of IR-fixed points in the model discussed: 

g , = Ae, 77* = AS, Vc„ A = {A7T) d/2 T{d/2)d{d - l)" 1 , (22) 
which exists in the physical region {g,r],c} > 0. Along this line 

7u = e, 72, = 25/3, (23) 

with no corrections of order e 2 , S 2 and so on. 

There are no new critical exponents in the model considered since the critical dimen- 
sionalities of (p, c, and v are just equal to the Kolmogorov's value —1/3 as it defined from 
the phenomenological assumptions (|j) and (^j), and the dimensionality of the field ip is 
equal to its canonical dimension since it has no anomalous corrections (7^ = 0), 

A v = A c = A„ = -1/3, A^ = d/2, A {i , H} = d. (24) 

The results expressed in the latter relation (^) can be applied to the study of relative 
diffusion, i.e., spreading of an admixture cloud consisting of a large number of particles, 
I6| . If we could label one of the particles in the cloud at time t = in the source of a 



frame x = 0, then the effective radius of the cloud R in the moment t > satisfies the 
relation 

R 2 = J dx x 2 (^ t (x,t)V(0,0)). (25) 

Taking into account that Ar = — 1, from (^) and (^) one can obtain that A [dR 2 /dt] = 
—2 — At, where A t is the Kolmogorov's dimensionality of time At = —2 + 71 = —2/3. 
It leads to dR 2 /dt ~ i? 4 / 3 , which was to be shown in the Richardson's Four-Thirds Law, 

From the first relation in (f24j) one can easily derive the well-known Kolmogorov's Law 



of Five-Thirds, 1TB 



We should stress that the indices ( ^4[) as well as their consequences (Richardson's and 
Kolmogorov's Laws) are the universal features of the model (|14"D, i.e., they do not depend 
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on the certain positions of fixed points of the RG-transformation of variables in the domain 
of IR-stability. Physically, this fact means that these properties do not depend on the 
parameters of the model {g, 77, c, u}. However, the scaling function T for the propagator G 
(and consequently the amplitude factors in the Kolmogorov's and Richardson's Laws, i.e., 
the Kolmogorov constant and the effective diffusivity coefficient) do depend on {g, 77, c, v}. 
We discuss it in more details later on. 

Another important note should be made on the role of parameter c for the static [i.e., 
of equal time) spectra. As we have mentioned before the parameter c determines the 
diffusion rate in the microscopic scale by means of porosity fluctuations. However, the 
results (|22| ) and (|24| ) demonstrate that c does not affect the universal properties of the 
static quantities of model but c can only change the amplitude factors of empirical laws. 

3.4 Solutions of RG-Equations, Asymptotics for the correlation 
function G and the Extinction Coefficient x 



Now we derive the solutions of RG differential equations (|T^) which give the IR - asymp- 
totics (r ^> w) for G. In renormalized variables the propagator G can be expressed in the 
form 

G- x = -k d R(s,g,ri,c,y,z), (26) 
where R is some function of dimensionless arguments 

s = k/M } y = uj/uM 2 } z = X 2 /M 2 } (27) 

which meets the RG-equation 

[-D s + f3 g d g + (3 v d v + (3 c d c - (2 - 7l ) A, + (2 7l - l2 )D z \ R = 0. (28) 



The equation ( |2lf ) can be solved and the function R can be found out. 

The general solution of (^) is an arbitrary function of the first integrals e = {g, f], c, 



y, z} which number is one less then the number of arguments of R in (|26|). They can be 
founded from the system of equations 

ds dg df] dc dy dz 



s (3 g (g) /3 v (g,v) Pc(g,v) 2/(2-71) 2(271-72)' 

supplied by some normalization conditions for e. We use the standard ones, 



(29) 



ei(s = l,ei) = Ci. (30) 



Traditionally, the first integrals associated with coupling constants are called invariant or 
running charges. In spite of the model (|14"D is of three charges the relations between them 
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are so simple that one can (at least up to one-loop order approximation) solve ( p9| ) 

[3 dx n dx 
lUS = J 9 JJx) = h JJx)' (31) 

These integrations can be brought about explicitly. For invariant charges it, in particular, 
leads to the expressions 

SM = es^ + Ml-^y ^ V) = 25^ + 3^(1-^) - (32) 

Taking into account the normalization condition Q3"0| ) and the relations (^Tj), for any 
solution R(s, g, i], c, y, z) we have 

fs e gY 2/3 fs 2S fjY 2/3 /c\- 2 
R(s,g,r),c,y,z) = I —J \~^~) \c) R ^ l ^9^^c,y,z). (33) 

Note, that the scaling function R(l, g, fj, c, y, z) is not fixed by RG-equations and is usually 
calculated in the framework of diagram technique. 

Under the renormalization group transformations the canonical degrees of freedom are 
replaced by the scaling ones including anomalies. For the scaling asymptotics s — > of 
(|33|) it is expressed in the fixation of the values of invariant charges on their values in the 
fixed points of RG-transformation ( p^) , 

9~(s,g)^g*, fj(s,g,r}) ->• 77*, s = k/M->0. (34) 

In the practical problem it would rather be convenient to express the asymptotics of cor- 
relation function G by virtue of bare parameters eo, using the relations of renormalization 
constants ( |T6| ) and fllTI) . Finally, one has the asymptotics s — > in the following form 

G- 1 ~ (35) 

where A = - (g) ~ 2/3 (j) ~ 2/3 (a) ~ 2 i?(l, r], c, y, 5). 

The essential feature of the model (|14D is traced in the dividing out of renormaliza- 
tion constants in the IR-asymptotics of Green's function as well as disappearing of the 
renormalization mass parameter M, which functions in ( |3"5] ) are commended to the bare 
parameters go,^,^,^. 

For the extinction index x from (p9|) we have 

-2 2 f S A 27l(g) 2 r , / 72(^) ,, R x 

X = X ex P / dx x ex P / dx . (3oJ 



In the analogous manner one can obtain the relevant scaling asymptotics: 

, , , \ y2 

,9* J \VoJ V c * 



Renormalization Group Approach 



16 



which corresponds to the following damping time-spectrum for the dynamic solutions G, 

G(t-t') ~exrj-a k? /3 \t-t'\, (38) 

where a = (j^j (^) (^) . Using (^), one can compute the asymptotic value of static 
concentration spectrum (with no time dependence) in the case of strong axial porosity 
fluctuations in the three dimensional space 

u(r) ~ -. (39) 
r 

This asymptotics replaces the gaussian one for the region much large then w, so that the 
ordinary gaussian spectrum for this region is to be replaced by k~ 2 for r ^> w, see Fig. 2. 

4 Renormalization and Critical Dimensions of Com- 
posite Operators c~ n (p n (x.) 



We have considered the model fll^f ) which corresponds to an unphysical version of the 
packed beds problem. Now we are going to generalize the developed technique to the case 
of nongaussian fluctuative component of the porosity field e. 

In the Section 2 we have explained that all amendments due to deviations from the 
gaussian distribution for e can be, in principle, taken into account (at least they can be 
well-estimated) as a series of power-like composite operators (c~ n (p n (x)) in the framework 
of renormalizat ion-group approach for the model (14). Now we calculate the precise 
critical indices for each power-like operator. 

The canonical dimensions of operators (c~ n <^ n (x)) are trivial, dp = —nd c + nd^ = 
0, Vn. From the Tab. 1. and ( p~3| ) one can see that the only diagrams with insertions of 
power-like operators which contain the superficial divergences have an arbitrary number 
of external "tails" of the c -1 <£>-type. However, at least one of such an external "tail" 
is attached either to the vertex %l)^(\d)ip or ip^(c~ l (pdt)ijj, so that at least one derivative 
d or dt appears as an extra factor in the diagram, and, consequently, the real index of 
divergence is necessarily negative. 

This means that all operators discussed require no counterterms, i.e., they are UV- 
finite. The same result can be readily reproduced by consideration of the Ward identities 
which express the invariance of the model ([14]) with respect to the field transformations 
(0). Generally speaking, the model (|T3) is invariant under the transformations flT2"|) but 
the operators (c _ri <^ n (x)) are not, so that they cannot have counterterms and cannot be 
UV-divergent. Finally, for the critical indices of power-like operators (</? n (x)) one obtains 
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the simple recurrent relation 

A[if n ] = -n/3. (40) 



Note that this relation was not initially clear, and it is a specific feature of the model (|14]) 
and those of similar type, [I3|. Analogously, one has A[c~ n (p n (x)] = 0, Vn. 

The latter result exhibits that no one term of power series for the fluctuative component 
of porosity field e can be omitted from the consideration in the equation (H) since their 
contributions are of equal importance. Thus, we need to consider the theory with all 
terms of the power series included in the action functional. 



5 On the Solutions for Arbitrary e 

The situation similar to that of the problem considered has place in the problem of 
chemically active scalar admixture [[OJ and passive scalar admixture convection |S{J . Each 
term of the power-like series happens to be significant in the region we are interested in, 
so that we cannot limit our consideration by neither the first term ~ c~ l (p as we did in 
the Sec. 3 nor any finite number of terms in the power series for a porosity field function 
(for example, (0)). 

Generally speaking, it would lead to the investigation of the problem of an infinite 
number of dimensionless coupling constants rj n , n — > oo, i.e., to the renormalized action 
functional (which is analogous to (0)) of the form 

Sr{^, v, V) = / dxdt [-|^D~V - fvD^ 1 v + (d t ij + Z 1 uAij + Z 2 vc~ 2 dfy+ 

+ c~ x <pd t ij) + En=2 c- n r] nV n d t ^/n\ + (vd)ip)] + C.C. 

(41) 

However, in contrast with the problems mentioned any of these new terms, c~ n r} n (p n d t i/j/n\, 
does not contain any UV-divergences in accordance with the general dimensional analysis 
of Sec. 3. 

This means that all the renormalization analysis of the theory (|14D is still also valid 
for fl4"ip . In fact, the coupling constants r] n are not the charges in theory (|4T|). Any of 
the critical indices of the theory flUp does not change in fl4"T| ) and no one new critical 
asymptotics appeared. 

The basic results (|35|), (|37|) , (p8|), ( [40"|) are still secure for any power series c~ n (p n in 
(p|) i.e., they are valid for any porosity function distributed in the microscopical range. 

The only difference between the various models of porosity for the certain critical 
asymptotics has place in the amplitude factors of Green's functions. They do depend 
from the certain initial values of the physical parameters of the problem. We believe 
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that it is the cause of a large dispersion of the various experimental measurements in the 
problem of the flow through the irregular packed beds. 

6 Conclusion 

The final conclusion is that the RG method applied to the problem of packed beds with 
the diffusion coefficient which can arbitrary depend from the concentration of marker 
in fluid flow demonstrates the existence of a scaling regime in the microscopical range 
where porosity is subject to the strong fluctuations. The relevant critical indices meet 
the empirical laws of Kolmogorov and Richardson for any porosity values. However, 
the nonuniversal quantities of the theory, i.e., the amplitudes of empirical laws are 
depended on the certain initial values of the parameters as well as the porosity fluctuations. 
The fluctuative porosity forms the damping time-spectrum for the dynamic correlation 
functions. Finally, the solution of nonlinear equation of diffusion for the concentration 
field w(x, t) is to be a traveling wave with an amplitude damping in space and time. This 
fact distinguishes the theory considered from the case of ordinary nonlinear diffusion as 
well as the case of turbulent mixing of passive advecture. 

We justified the result of |4]] on the inclusion of a "time- diffusive" term in the equation 
(0) for accounting of porosity fluctuation and justified it for an arbitrary stochastic process 
originating the porosity fluctuations. 
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Table 1. Canonical Dimensions of Quantities of the Theory 



F 


v,c,ip 




9o 


Vo 


9,V 


m, fi, M 




d% 


-1 


d/2 


e 


25 





1 


-2 


d F 


1 

















1 


dp 


1 


d/2 


e 


25 





1 






Renormalization Group Approach 



21 



CAPTION FOR FIGURES 
FIGURE 1. On the definition of porosity. 

The porosity is directly related to the size of the pores relative to the matrix. If P be a 
point inside a porous medium on the column axis surrounded by a layer of volume V(Zj), one 
can define a ratio, Si = ™ wn i cn Vv(k) is the volume of void space within volume V(U). 

FIGURE 2. The deviation of dimensionless concentration spectrum from gaussian for 
r w. 

Strongly developed fluctuations of the porosity and velocity fields lead to anoumalously small 
fading off of concentration on the distances r 2> w larger than the size of a pore. 



